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THE  ROLE  OF  FINITE  PARALLEL  LENGTH 
ON  THE  STABILITY  OF  BARIUM  CLOUDS 

I .  INTRODUCTION 

In  the  disturbed  ionosphere,  Rayleigh-Taylor  [Scannapieco  and  Ossakow, 
1976;  Ott,  1978;  Ossakow  et  al.,  1979;  Zalesak  and  Ossakow,  19801  and 
gradient  drift  instabilities  [Simon,  1963;  Linson  and  Workman,  1970; 
McDonald  et  al.,  1980,  1981]  are  considered  to  be  primary  sources  for  fluid 
structuring.  These  instabilities  are  predicted  to  be  aligned  along 
magnetic  field  lines  with  the  result  that  analysis  and  numerical  simulation 
has  emphasized  the  flute  approximation,  which  neglects  the  explicit 
dependence  of  modes  on  the  coordinate  parallel  to  the  ambient  magnetic 
field;  when  considered,  it  has  been  taken  into  account,  in  a  gross  sense, 
by  averaging  plasma  parameters  over  magnetic  field  lines,  but  without 
considering  the  implications  on  the  ionosphere  of  accompanying  mode 
variation  along  magnetic  field  lines.  One  exception  to  this  generality  is 
the  paper  by  Goldman  et  al.  [1976],  which  calculates  eigenmodes  in  the 
electrostatic  approximation,  recognizing  that  modes  must  vary  along  the 
magnetic  field  lines  as  one  moves  away  from  the  source  of  instability,  an 
artificial  plasma  cloud. 

Artificial  plasma  clouds  in  the  ionosphere  have  a  finite  spatial 

extent  that  can  influence  the  development  and  properties  of  plasma 

Instabilities.  In  this  paper,  particular  emphasis  is  placed  on  determining 

certain  effects  on  the  gradient  drift  Instability  resulting  from  the  finite 

length  of  plasma  clouds  along  the  geomagnetic  field.  It  is  demonstrated, 

with  a  rather  simple  plasma  geometry,  that  finite  plasma  length  implies 

parallel  currents  and  electric  fields  which  contribute  to  the  development 

of  eigenmode  structure  along  the  geomagnetic  field.  Specifically ,  the 

finite  field-line  integrated  electron  density  of  ionospheric  plasma  clouds 

can  play  a  role  In  reducing  the  growth  rate  of  striation  instabilities. 
Manuscript  approved  June  27,  1984.  1 


This  point  Is  subsequently  quantified.  The  finite  parallel  length  of 
ionospheric  plasma  clouds  tends  to  favor  the  growth  of  strlatlons  with 
short  perpendicular  wavelengths. 

In  Section  II,  the  calculation  model  and  general  equations  are 
discussed.  Section  III  contains  a  derivative  of  the  dispersion  equation. 
Section  IV  describes  the  quantitative  evaluation  of  the  model  for 
parameters  appropriate  to  barium  releases.  Concluding  remarks  are  given  in 
Section  V. 


II.  GENERAL  EQUATIONS 

We  first  derive  a  set  of  nonlinear  equations  to  describe  the  evolution 

A 

of  a  cold  plasma  cloud  (T  “0)  in  a  uniform  magnetic  field  B  ■  z  with  a 

A 

background  neutral  wind  Yn  ”  Vn  X  fsee  *al*  For  we 

consider  only  low  frequency  3/3t  «  motion  of  the  cloud  and  take  the 

electron  collisions  to  be  sufficiently  weak  so  that  v  /Q  «  1  hut 

e  e 

allow  to  be  arbitrary.  The  collision  frequency  and  gyrofrequency  of 

the  a  species  are  given  by  and  respectively.  In  this  limit  the 
fundamental  equations  of  our  analysis  are  continuity,  momentum  transfer, 
charge  neutrality  and  Ampere's  law: 


3n 

— — -  +  7  <  [nv  1*0 

3t  ^  ~or 
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0  *  -  eE  -  —  v  x 
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where  the  variables  have  their  usual  meaning.  We  take  the  electric  and 
magnetic  fields  to  be  represented  by  potentials  as 


1  3Az  " 

i  "  -  -7^  2 


B  *  B-  z  +  7A  x  z 
z 


where  <£  is  the  electrostatic  potential  and  ^  is  the  vector  potential, 
associated  with  the  magnetic  field  produced  by  the  self-consistent  plasma 
currents.  We  consider  only  Az  since  and  we  assume 


I 7A  x  z|  «  B_. 
z  u 


The  electron  cross-field  motion  is  given  by 


Sei  ‘  '  B  V  *  2 


while  the  parallel  motion  is  given  by 


e  »  1 

v  ■  -  I  f  b  •  7  + - 1 

ell  ra  v  1 '  v  c  3 1 

e  e 


where  v  *  v  .  +  v  and  b  *  B/B„. 
e  ei  en  -  0 
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The  ion  cross-field  motion  is  given  by 


the  charge  conservation  equation,  and  (14)  is  Ampere's  law.  Thus,  (12)- 
(14)  provide  a  complete  description  of  the  evolution  of  a  three- 
dimensional,  cold  plasma  cloud. 

We  will  only  consider  the  linear  stability  of  a  two-dimensional  cloud 
which  is  localized  both  along  and  across  the  magnetic  field  Bq:  nc  « 

nc(x,z)  with  |x|  <  Xq  and  |z|  <  Zq.  The  background  plasma  is  taken  to  be 
uniform  throughout  all  space.  The  equations  describing  this  two- 
dimensional  equilibrium  are  given  by 


3  A  * 

3t  e  3z  n  c  3 1  3z 


2 

Ac  3_  3 d>  .  c  3  3  Az  _ 

W.  B  3x  n  3x  4ire  3z  2 
i  3x 


32A 


3x 


a  a  * 

4tt  .  1  _ z  34 , 

cn  ’c  3t  3z-' 

e 


(15) 

(16) 

(17) 


From  (15)-(17)  we  find  that  the  equilibrium  is  given  by  A  *4*0  with 

z 

nc(x,z)  an  arbitrary  function.  For  simplicity  we  consider  the  plasma 
density  to  be  given  by  [similar  to  that  used  by  Sperling  (1983a)] 


n(x,z) 


n<  -  nc(x)  +  nb 


1*1 

1*1 


<  z, 


>  z. 


(18) 


where  the  subscripts  c  and  b  refer  to  cloud  and  background,  respectively, 
[see  Fig.  lb]. 

We  note  that  if  we  had  retained  a  finite  plasma  temperature  then  the 
density  equation  would  be 


3n  _  d—  T  3n 

3t  3z  m.v.  3z 

i  in 


(19) 
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where  T  *  Tg  +  and  we  have  taken  3/3x  -  0.  Equation  (19)  describes  the 

diffusion  of  the  cloud  along  z  with  a  diffusion  coefficient,  Dz  =■  T/m^v^. 

Thus,  in  general  (15)-(17)  will  not  have  equilibrium  solutions 

when  3/3z  *  0  and  T  *  0.  However,  in  calculating  the  growth  rates  of 

unstable  modes  with  growth  times  which  are  short  compared  with  the 

2 

diffusion  time  tz  ■  Zq/Dz>  we  would  expect  the  evolution  of  the  equilibrium 
to  have  very  little  influence  on  the  stability  calculation  for  this 
situation.  Consideration  of  finite  temperature  effects  will  be  deferred  to 
a  later  report. 


III.  LINEARIZED  EQUATIONS  AND  DISPERSION  RELATION 

To  find  the  influence  of  the  parallel  dynamics  on  the  instability,  we 

linearize  (7)-(9)  where  the  perturbed  quantities  vary  as  p  -  p(z)  exp  [yt  + 

ik^yj.  After  eliminating  the  equation  for  n  algebraically,  we  obtain  two 

coupled  differential  equations  for  A  and  |, 

z 


(20) 


,  3A 
_  1  z 

-  yoD  —  — — 
r  c  3z 


(21) 


where  y_  =  -  n;v  /n  5,  n'  -  3n./3x,  y  -  y  +  ik  V  v.  /Cl.,  D  =v  c2/w  2  is  the 
-j  snss  s  yn  m  t  r  e  pe 

resistive  diffusion  coefficient,  and  i  •  ~l  Q,/v  v,  5. 

e  i  s  in 

Prior  to  solving  (20)  and  (21)  for  the  density  profile  given  by  (18), 
we  first  consider  a  cloud  of  infinite  extent  [z^  *  »  in  (18) J  and  Fourier 
expand  modes  parallel  to  30 ,  i.e.,  p(z)  -  p  exp  [ik^zj.  We  consider  the 
local  approximation  [k^n'/nj  ^  >>  lj  so  that  (20)  and  (21)  can  be  solved 


algebraically.  This  allows  comparison  with  previous  results  [Sperling, 
1983b;  Basu  and  Coppi,  1983]  and  insight  into  the  effects  of  finite 
parallel  wavelengths.  The  local  dispersion  equation  is  given  by 


(y  -  Y0J(y  +  kjDrJ 


-  yokzDr 


We  note  that  the  RHS  of  (22)  can  be  expressed  as  -  yk  5,  where  V.  ■ 

z  A  in  i  a 

1/2  - 

B/^Aimm^J  is  the  Alfven  velocity.  This  form  explicitly  shows  the 

* 

coupling  to  an  Alfven  wave  which  is  the  dominant  finite  k2  effect.  For 

simplicity  we  assume  y  »  k  V  ( v.  /n,  )  and  solve  (22)  in  two  limits:  y  » 

y  n1,  in  l '  1 

2  2 
kyDr  (electromagnetic  limit)  and  y  «  (electrostatic  limit).  It  is 

found  that 


1  -  - 1  - 


It  »  Vt) 


-1 

(1  +  kjAiyJ 


[y  «  k  u  ) 
v  '  y  r 


where  y  *  y/yn,  k  ■  k  L  ,  k  ■  k  z  ,  and  L  and  z  are  the  perpendicular 
^  y  y  i  z  z  r  ^  * 

2  2 

and  parallel  resistive  length  scales  given  by  -  D^/y^  and  z^  »  aD^/y^. 

In  the  ideal  limit  (Dr  *  0),  (23)  indicates  that  electromagnetic  effects 

*2  2  2 

stabilize  the  instability  when  Y.^  >  1  (i.e.,  kzVA  >  y^Vj  -5^)  by  coupling  to 
* 

an  Alfven  wave.  This  result  is  similar  to  that  of  Basu  and  Coppi  [1983] 
who  investigated  the  collisional  Rayleigh-Tavlor  instability.  On  the  other 
hand,  (24)  indicates  that  in  the  electrostatic  limit  there  is  only  a 
reduction  in  growth  rate.  This  is  because  resistive  diffusion  across  the 


7 


magnetic  field  lines  is  sufficiently  rapid  to  dissipate  the  Alfven  wave. 

Finally,  in  the  limit  k2  +0  we  recover  from  (22)  the  standard  result 

Y  ■  Yq  ~  ik  vin/^i  J »  w^ere  the  real  frequency  is  caused  by  the 

equilibrium  electric  field  in  the  x  direction;  in  a  frame  of  reference 

moving  with  the  electrons,  however,  the  mode  has  no  real  frequency.  A 

2 

second  mode,  which  is  damped  is  given  by  y  *  -k^D^  ~  ik^V^  J  t^hu  et 

al.,  1978;  Sperling,  1983b]. 

For  a  general  profile  n(x,z)  the  coupled  equations  for  $  and  A^  must 

be  solved  subject  to  the  boundary  conditions  A,  A  +  0  as  I z  I  -*■  ®.  For  the 

z 

step  profile  for  n(x,z)  given  in  (15),  the  solutions  to  (20)  and  (21)  in 
the  region  |z|  >  Zq  can  be  written  as  plane  waves, 


Az  =■  Az>  exp  [-  k>  I  z  |  J ,  (25a) 

}  *  exp[-  k> | 2 1 J,  (25b) 


with 


where 


+  Vr> 


JA 


z> 


k>c 


(I  i_] 

v  y  aRO 


k2D 


y  r 


Y  +  k'D 


(25c) 


( 25d ) 


R 


(25e) 


and  the  subscript  >  on  a  given  parameter  indicates  that  it  is  to  be 


evaluated  in  the  region  |z|  >  z^.  Mote  that  the  solutions  which  diverge 

as  z  +  +  «  have  been  omitted  in  (25a)  and  (25b).  The  parameter  R  is  a 

measure  of  the  electrostatic  or  electromagnetic  nature  of  the  mode.  For 
2 

kyD^  »  y,  we  note  that  R  -  1  and  the  mode  is  essentially  electrostatic. 

2 

In  the  opposite  limit  k^D^  «  Yj  the  mode  is  electromagnetic  and  R  «  1. 

In  the  region  |z|  <  z^ ,  the  solutions  for  $  and  Az  are 


K  "  K<  sltHk<zJ  '  K<  C08ik<zJ 


(26a) 


$  *  cosjk^zj  +  sin^k^zj 


(26b) 


with 


1,2 

z< 


k<c 


(26c) 


(26d) 


To  complete  the  dispersion  relation,  we  must  match  the  various  plane  wave 
solutions  at  z  *  +  Zq  .  The  appropriate  matching  conditions  are  obtained 
from  (20)  and  (21).  Integrating  these  two  equations  across  the 
discontinuity  in  the  density  at  z  =  +  zn,  we  find  that  $  and  Az  must  be 
continuous.  For  the  even  $  solution  *  0),  we  obtain  the  dispersion 

relation  for  the  growth  rate  y, 


where  m  is  an  integer.  The  dispersion  relation  for  the  odd  <J>  mode 

( =  0 )  is  similar  to  (27)  except  tan-^  is  replaced  by  -  cot 

In  general,  the  dispersion  equation  (27)  has  an  infinite  number  of 

solutions  for  a  given  set  of  physical  parameters,  corresponding  to 

eigenmodes  with  an  increasing  number  of  modes  (m)  along  z.  The  dispersion 

equation  (27)  can  be  solved  numerically  for  arbitrary  values  of  the 

background  and  cloud  parameters.  However,  to  gain  an  understanding  of  the 

general  scaling  of  the  growth  rate  y  with  the  parallel  extent  of  the  cloud, 

it  is  useful  to  solve  (27)  analytically.  To  do  this  we  make  a  number  of 

simplifying  assumptions.  We  consider  only  the  lowest  order  mode  (i.e. ,  m  * 

0  which  implies  0  <  k^zQ  <  ir/2;  it  is  easily  shown  that  this  mode  has  the 

largest  growth  rate);  take  v.  /n.  to  be  small  so  that  5=1  and  the  real 

in  i 

frequency,  k^V^v^/fl^,  can  be  neglected;  and  finally,  take  all  parameters 
but  the  density  to  be  the  same  inside  and  outside  of  the  cloud.  In  this 
limit  the  expressions  for  k^  reduce  to 


(28a) 


(28b) 


We  separate  our  analysis  into  two  separate  cases  v  .  >  -j  .  In  the  limit 

ei  <  en 

v  .  >>  v  ,  the  resistive  diffusion  coefficient  is  continuous  across  the 
ei  en  r 


cloud  boundary  and  Rv  =  Rx. 


K  > 


With  these  assumptions  the  dispersion  equation  (27)  becomes 


k<z0  *  tan_1U— V  ~)l/2\  (29a) 

1  -  Y  < 

with 

*<  4  ‘  It1  '  +  y^O2.  (2,b) 

~  2  2  2 
and  z-  ■  z-/z  >  k  *  k  L  where  L  ■  D  v/y  and  z  *  •  The  dispersion 

0  Or’y  yr  r  r>'o  r^>r  K 

relation  in  (29)  is  now  a  function  of  only  three  parameters:  zQ,  n^/n^  an*1 

y* 

In  the  limit 


[(1  -  Yj/Yj(n</n>J  «  1, 


(30) 


the  arctan  function  approaches  ir/2  so  that  k^z^  =  ir/2  is  the  dispersion 
equation.  Furthermore,  since  n^/n^  *  (n£  +  n^J/n^  >  1,  the  inequality  in 
(30)  can  only  be  satisfied  for  y  »  1.  Thus,  the  growth  rate  is  given  by 


and 


-  .  vi_n> 

Y  "  1  *2  n 

4z0  < 


for 


y  »  k2 

y 


*  n  .  *2  n>  1  ,-l 
Y  Vi  +  —  T.tz) 

4z0  <  ky 


for 


y  «  k2 

y 


(31a) 


(31b) 


M 


» 
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We  note  that  (31a)  and  (31b)  can  be  compared  to  the  local  growth  rates  (23) 
and  (24)  If  we  define  an  effective  k 


“2  .  n2  n> 

Czeff  =  72  n  .  ’ 

< 


“2  “ 2 

Thus,  taking  kz  -  kgeff  in  (23)  and  (24)  we  recover  (31,a)  and  (31b).  For 
m  ^  0  modes,  becomes  (m  +  1/2)  it  Mzg^n^/n^J  indicating  a  spectrum  of 
modes  parallel  to  B.  Finally,  since  y  ~  1  from  (31a)  and  (31b)  we  can 
write  a  generalized  dispersion  relation  given  by 


■y  - 1  - 

,  2  y  >  n . 


(33a) 


and  is  valid  for  [based  on  (30)] 


I!  +  kyjZo  »  I* 


(33b) 


In  the  opposite  limit,  i.e., 


U1  *  Y  J /  Y  J  (n</n>  J  »  1, 


arctan  function  approaches  k^z^.  The  dispersion  relation  becomes 


l1  ”  YJU  +  ky2J1/2  "  Y  n>/n<zQ< 


The  solution  of  the  equation  is 


Y  *  1  -  '.n>/n< JZq1  [1  +  \2)'U2  1 


(36a) 


f 


L 


¥ 


i 

b. 
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[n>/n<j2  «  (1  +  ky2Jz02  «  1,  (36b) 

^  A 

where  the  inequalities  follow  from  (34)  and  the  condition  y  »  1.  When  y  « 
1,  (35)  becomes  a  quadratic  equation  with  the  solution 

y  -  z02(n</n> J2[l  +  [1  +  4  ky2z'2n2/n2 J 1/2 j /2,  (37a) 

which  is  valid  for  (from  Y  «  1) 

l1  +  \2fz02  «  n>/"<  •  (37b) 

For  zQ  very  small  the  growth  rate  approaches  zero  as 

Y  s  kyzon</n>*  (38) 

A 

Equation  (38)  also  implies  that  the  growth  rate  increases  with  ky  in  this 
region.  The  physical  reason  for  this  behavior  will  be  discussed  shortly. 

The  growth  rates  and  inequalities  in  (33),  (36)  and  (37)  can  be 
summarized  rather  succinctly  in  the  (n^/n^J  -  z^  phase  space  plot  shown  in 
Fig.  2.  Note  that  n^/n^  *  n^/(n^  +  ncJ  <  1*  Regions  I,  IT  and  III 
indicate  the  range  of  validity  of  the  growth  rates  in  (33),  (36)  and  (37), 
respectively.  F.xpressions  for  the  growth  rate  in  these  three  regions  are 
given  in  Table  I.  For  large  z^ ,  in  Regions  I  and  II,  the  finite  length  of 
the  cloud  along  3  has  very  little  influence  on  the  growth  rate 
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1. 


and  y  -  Y/ Yq  -  1-  For  small  ,  In  Region  111,  y  <<  1  so  that  the  finite 
extent  of  the  cloud  strongly  reduces  the  growth  rate.  As  ky  Increases, 
Region  III  shrinks  In  size  so  that  the  growth  rates  for  short  wavelength 
modes  are  less  affected  by  the  parallel  dynamics  than  are  the  growth  rates 
for  long  wavelength  modes.  Nevertheless,  (25d)  shows  that  the  short 
wavelength  modes  are  more  localized  along  the  magnetic  field  than  long 
wavelength  modes  [Sperling  and  Glassman,  1983;  Sperling,  1984]. 

A 

In  Fig.  3  we  schematically  show  the  growth  rate  y  as  a  function 

A  A 

of  Zq  for  n^/n^  and  ky  «  l  held  fixed.  In  Region  III  the  growth  rate 
first  increases  linearly  with  z^  and  then  quadratically  until  y  -  1>  where 
it  enters  Regions  II.  The  transition  from  y  ~  to  y  -  zq  in  Region  III 
is  a  consequence  of  the  change  in  character  of  the  mode  from  being 
dominantly  electrostatic  (y  <  k^D^J  f°r  zq  ^  'tyn>^n<  to  electromagnetic 
U  >  k^DrJ  for  zQ  >  kyn>/n<. 

A 

The  dependence  of  the  growth  rate  on  ky  can  also  be  readily  obtained 

A 

from  Fig.  2  and  Table  I.  For  parameters  z^  and  n^/n^  such  that  modes  with 

A  AAA 

ky  «  1  are  in  Regions  I  and  II  (z^  >  n^/n^J,  y  -  1  for  nil  ky  since 

A 

increasing  ky  simply  pushes  the  mode  further  into  Regions  I  and  II.  For 
the  case  where  modes  with  ky  <<  1  fall  in  Region  III  (z^  <  n^/n^J,  the 
dependence  of  y  on  ky  is  more  interesting.  In  Fig.  4  the  growth  rate  is 
shown  versus  k  with  nN/nx  and  z.  held  fixed.  ror  small  k  the  mode  falls 
in  Region  III  and  has  a  growth  rate  nearly  independent  of  ky.  When  ky  > 
z^n^/n^,  the  mode  becomes  electrostatic  and  the  growth  rate  Increases 
with  k  until  k  -  nv/n.z-  when  y  '  1  and  the  mode  enters  Region  II.  We 
again  conclude  from  this  figure  that  the  growth  rates  for  short  wavelength 


modes  are  less  affected  by  the  parallel  dynamics  than  the  growth  rates  for 


long  wavelength  modes. 
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Up  to  this  point  we  have  discussed  only  the  case  where  »  vgn. 

The  dispersion  equation  (27)  can  be  solved  in  a  similar  manner  in  the 

opposite  limit  «  ven-  In  this  limit  the  resistive  diffusion 

coefficient  is  not  continuous  across  the  cloud  boundaries  at  +  Zq.  The 

results  differ  only  slightly  from  those  just  presented  so  we  skip  the 

details  and  simply  present  the  analogues  of  Fig.  2  and  Table  I  for  this 

case.  In  Fig.  5  we  show  the  [n^/n^J  -  Zq  phase  space  plot  showing  the 

three  regions  of  the  instability  for  vgn  »  \>e^.  The  growth  rates  for 

these  three  regions  are  listed  in  Table  II.  The  only  difference  between 

phase  space  plots  in  Figs.  2  and  5  is  the  boundary  between  Regions  I  and 

II ,  which  falls  at  larger  zrt  when  v  »  v  . .  The  growth  rates 

')  en  el 

for  v  v  .  differ  only  in  Region  I.  Since  y  -  1  in  both  Regions  I  and 

en  <£  ei  ® 

II,  these  differences  are  not  particularly  significant  so  the  previous 

discussion  of  the  instability  in  the  limit  v  ,  »  v  also  applies  to  the 

ei  en 

opposite  limit. 

We  have  shown  that  the  finite  length  Zq  of  the  plasma  cloud  reduces 
the  growth  rate  of  the  gradient  drift  instability  compared  with  its  value 
when  Zq  is  infinite.  The  growth  rates  of  long  wavelength  modes  are  more 
strongly  reduced  than  those  of  short  wavelength  modes.  The  essential 
physics  which  underlies  these  results  can  be  understood  by  integrating  (21) 
for  the  perturbed  potential  i  along  z, 

»  z0 

Y  /  dz  n(z)$(z)  -  yQ  J  dz  n(z)<i>(z)  =  2Y0n<zf)$<,>  (39) 

"Z0 

where  we  have  again  taken  small  and  we  have  assumed  that  >(z)  = 

const  in  the  region  ] z I  <  z^ ,  which  is  approximately  correct  for  the  lowest 


I 


I 


order  mode. 


The  integral  on  the  left  side  of  (39)  represents  the 


integrated  Pederson  conductivity  along  the  entire  field  line,  weighted  by 


the  potential  The  integral  on  the  right  side  of  this  equation  is  the 
integrated  Pederson  conductivity  of  the  cloud.  Carrying  out  the  remaining 


integral,  we  find 


Y  “  Y0n<z0 Vtn<z0*<  +  n>k> 


*  Yon<Vtn<zo  +  n>k>  )' 


where  we  have  used  the  results  of  our  previous  calculation:  ^  and 

^  >  0  when  k^|z|  >  1.  Substituting  the  expression  for  k>  in  (25d)  into 


(40),  we  obtain  the  dispersion  relation 


Y  -  Yn 


.-1,  1/2 

n<zQ  +  n>ky 


When  the  result  is  written  in  the  normalized  variables  presented 
previously,  it  reduces  precisely  to  the  dispersion  relation  presented  in 
(35).  The  reduction  of  the  growth  rates  due  to  finite  Zq  occur  because  the 
integrated  Pederson  conductivity  over  the  extent  of  the  mode  along  z 
becomes  comparable  to  that  of  the  cloud.  The  distance  the  mode  extends 
along  z  k^  J  increases  as  k^  decreases  [see  (25d)]  so  that  long 
wavelength  modes  have  greater  reductions  In  their  growth  rates.  The 


electromagnetic  effects  prevent  the  mode  from  extending  an  Infinite 
distance  along  the  magnetic  field  line  as  k^  *  0.  Tn  this  limit  R^  ♦  0  and 
(41)  becomes 


Vo 


oD_ 


(42) 


(zo  +  n>^“rJ> 


so  that  the  growth  rate  Is  independent  of  ky  as  shown  In  Fig.  4. 


IV.  RESULTS 

The  plasma  model  outlined  In  Sections  II  and  III  Is  rather  simple,  and 

admittedly  does  not  include  all  the  intricacies  and  inhomogeneities  of  the 

real  ionosphere.  Nevertheless,  the  model  does  show  that  the  growth  rate  Is 

a  function  of  the  parallel  length.  We  now  present  quantitative  results  of 

our  model  for  typical  Ionospheric  parameters.  In  particular,  we  take 
7-3  5-3 

n^  *  10  era  ,  n^-  10  cm  ,  B  »  0.4  G,  T^  -  Tg  -  0.1  eV,  and  consider  two 

values  of  Yq '  Yq  *  0.01  sec  *  and  Yq  *  0.10  sec”*.  The  Ions  for  |z|  <  z ^ 

are  assumed  to  be  barium  (i.e.,  m^  »  137  m^,  with  nip  the  proton  mass)  and 

the  Ions  for  |z|  >  zn  are  assumed  to  be  air  with  m.v  ■  20  m  .  The  mass  of 

u  1?  p 

an  atmospheric  neutral  particle,  is  also  assumed  to  be  20  mp.  The 

9  -3 

neutral  density,  n^  ■  n^  »  8.3  x  10  cm  ,  is  appropriate  to  an  altitude 
of  ■»  200  km  [Knapp  and  Schwartz,  1975].  The  collision  frequencies  are 
determined  from  [Braginskii,  1965;  Xilb,  1977) 

*  2.5  x  lO^m^n^  (42a) 

vt>  -  2.5  x  10 13  m0n0  (42b) 
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_q 

v  ...  v  -  0.51  v  +  1.9  x  10  nA 

e<(>)  ei<(>)  0 


-3 

u  .....  *  1.1  x  10  n  ..... 
ei<(>)  e<(>) 


(42c) 


(42d) 


In  (42)  all  parameters  are  In  cgs  units.  Based  upon  these  parameters  and 
equations  we  can  calculate  the  perpendicular  and  parallel  resistive  scale 
lengths  which  are  given  In  Table  III  for  Yq  ■  0.01  sec  *  and  0.1  sec  *. 

An  important  result  to  be  considered  Is  the  boundary  between  Regions 
IT  and  III  In  Figs.  2  and  5.  The  boundary  is  defined  by  the  equation 


>  -21/2“ 
—  -  (1  +  k  V 

n^  y  0c 


and  marks  the  transition  between  strong  and  weak  growth.  For  z^  <  z^c  the 

theory  predicts  very  weak  growth,  and  it  is  clear  that  the  long  wavelength 

modes  (k^  «  1)  have  the  weakest  growth  for  a  given  value 

of  Zq  and  n^/n^.  In  Fig.  6  we  plot  Zqc  (km)  vs.  k^  (km  )  based  on  (43)  for 

n^/n^  «  0.01  and  the  parameters  described  above  and  in  Table  III.  The 

regions  below  each  curve  indicate  weak  growth,  while  regions  above  each 

curve  indicate  strong  growth.  As  expected,  the  value  of  z0c  is  largest  for 

ky  «  1.  Also,  we  note  that  the  dependence  of  zr)c  on  ky  is  stronger  for 

smaller  values  of  wavenumber  (k  <  l  km  S  and  . 

y  0 

In  Fig.  7  we  plot  Y  *  y/ Yq  vs.  ky  (km-^)  for  the  parameters  given 
above  and  different  values  of  z^  .  The  solid  lines  are  based  on  Yq  “  0.10 
sec-^  while  the  dashed  lines  are  for  Yq  ■  0.01  sec  ^ .  These  figures  are 
based  upon  (35).  Again,  these  curves  quantitatively  reflect  the  growth 
rate  dependence  on  ky  and  Zq  which  is  presented  qualitatively  in  Fig.  4. 
For  k  >  2  km  *,  the  normalized  growth  rate  y  is  independent  of  Zr, .  This 
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is  consistent  with  Fig.  6  which  shows  that  z^c  is  independent  of  y^  for 

k  >  2  km  ^ .  On  the  other  hand,  for  k  <  2  km  ^  the  difference  between  the 
y  v 

curves  becomes  larger,  especially  for  large.  Again,  this  is  consistent 
with  Fig.  6  which  indicates  a  strong  dependence  of  zqc  on  y^  for  ky  < 

2  km-^.  It  is  clear  from  Fig.  7  that  the  finite  parallel  extent  of  the 
cloud  favors  the  growth  of  short  wavelength  modes. 


V.  CONCLUDING  REMARKS 

We  have  presented  a  linear  analysis  of  the  gradient  drift  instability 
appropriate  to  ionospheric  plasma  clouds  of  finite  spatial  extent,  both 
perpendicular  and  parallel  to  the  ambient  magnetic  field.  Based  on  a 
simple  model  for  the  cloud  [see  (18)],  we  find  that  the  parallel  extent  of 
the  cloud  along  B  can  significantly  alter  the  growth  rate  of  the 
instability.  Specifically,  we  find  that  long  wavelength  modes  (ky  small) 
tend  to  have  much  smaller  growth  rates  than  short  wavelength  modes  (ky 
large).  This  reduction  in  growth  due  to  finite  ^  occurs  because  the 
integrated  Pedersen  conductivity  of  the  background  plasma  over  the  extent 
of  the  mode  along  z  becomes  greater  than  that  of  the  cloud  [see  (40)].  The 
distance  the  mode  extends  along  z  increases  as  k^  decreases  fsee  ( 2 5d ) ]  so 
that  the  long  wavelength  modes  are  more  strongly  affected.  A  more 

detailed  numerical  analysis  of  this  effect  is  given  in  Sperling  and 
Glassman  (1983).  We  mention  that  these  results  also  aoply  to  the  F  x  R 
gradient  drift  [McDonald  et  al.  ,  1980,  1981]  and  Rayleigh-Tavlor  [Ossakow 
et  al.,  1979;  Sperling,  1982]  instabilities. 

A  question  which  naturally  arises  is,  how  does  the  growth  rate  change 
as  the  cloud  diffuses  along  the  magnetic  field  or,  more  generally,  as  the 
cloud  evolves  in  time?  It  is  tempting  to  conclude  that  the  growth  rate 
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simply  increases  with  zg  as  the  cloud  diffuses  along  z  as  shown  in  Figs.  3 
and  6.  However,  as  the  cloud  diffuses  along  z  we  expect  the  electron 
density  inside  the  cloud  to  decrease  such  that  n^zg  -  constant.  Loosely 
speaking,  the  total  number  of  cloud  electrons  along  a  given  magnetic  field 
line  should  not  change,  or  at  least,  not  change  significantly.  Since  the 
growth  rates  of  the  instability  in  the  weak  (Region  III)  and  strong  (Region 
II)  growth  regions  depend  only  on  the  product  n<,zg,  parallel  diffusion 
along  z  does  not  cause  the  growth  rate  to  increase.  The  expansion  of  the 
cloud  along  z  also  does  not  enable  the  mode  to  change  from  weak  to  strong 
growth.  This  is  clearly  seen  from  (43)  which  indicates  that  the  boundary 
between  Regions  II  and  III  does  not  change  in  time  as  long  as  n^Zg  - 
constant.  Thus,  we  conclude  that  just  diffusion  of  the  cloud  along  z  does 
not  change  the  growth  rate  in  an  important  way. 

On  the  other  hand,  the  gross  evolution  of  the  cloud  involves  not  only 
diffusion  along  z  but,  for  example,  steepening  of  the  "backside"  of  the 
cloud  which  increases  the  growth  rate  Yg*  As  shown  in  Fig.  6,  for  a  given 
value  of  zg  more  modes  are  in  the  strong  growth  regime  when  Yg  is  larger. 
Thus,  the  parallel  dynamics  will  significantly  change  the  spectrum  of 
unstable  modes  as  the  cloud  evolves  in  time.  This  feature  may  be  related 
to  the  observation  of  the  delayed  onset  of  striations.  That  is,  the  onset 
of  striations  in  barium  clouds  occurs  several  minutes  after  the  detonation 
of  the  barium  release.  Specifically  for  larger  barium  clouds,  like  Spruce 
[Linson  and  Meltz,  1972],  the  onset  time  for  striations  is  estimated  to  be 
generally  greater  than  approximately  5  minutes  [McDonald  et  al.  ,  1980]. 

Linson  and  Meltz  (1972)  have  estimated  that  the  Spruce  and  Olive  ion  clouds 
had  an  extent  approximately  30-40  km  extent  along  the  magnetic  field  at 
about  5  minutes.  This  parallel  extent  corresponds  to  z0  ■  15-20  km.  It  is 


possible  then  that  plasma  clouds  need  to  evolve  in  time  so  that  they  have  a 


sufficiently  steepened  "backside",  as  well  as  a  sufficient  length  along  3, 
before  the  striations  can  occur  rapidly. 
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TABLE  I :  GROWTH  RATES  f  v  ,  »  v  1 
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TABLE  II:  GROWTH  RATES  fv  J  «  v  | 
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Region 

Growth  Rate  j 

— 

(1  +k  /rA 

Regions  of  validity  for  the  growth  rates  given  by  (33),  (36), 
and  (37).  The  growth  rates  are  listed  in  Table  I  and  are  based 


upon  the  assumption  that  v  »  v  . 


'ei 


»  v, 


en 


k  !!>  — 

v  n<  A  nc 


Schematic  of  the  growth  rate  y  vs.  parallel  length  for  fixed 

A 

n^/n^  and  k  .  This  curve  is  based  on  the  assumption 


that  v  .  >>  v 


<  ^ 


Regions  of  validity  for  the  growth  rates  in  the  regime 

v  .  <<  v  .  The  growth  rates  are  listed  in  Table  II.  This 
el  en 

curve  is  similar  to  Fig.  2  but  the  boundary  between  regions  I 
and  II  is  shifted  to  larger  zn. 


y0  =  0.01  sec  1 


Plot  of  "critical"  parallel  length  z  (km)  vs.  perpendicular 

Oc 


-1. 


wavenumber  k^  (km  )  based  upon  (43)  for  typical  ionospheric 


parameters.  The  parameters  used  are  given  in  the  text  (Section 
IV).  These  curves  denote  the  boundary  between  Region  II  [above 
curves  -  strong  growth  ( y  —  1 )  ]  and  Region  III  [below  curves  - 

A 

weak  growth  (y  <<  !)]• 
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UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOflETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01 CY  ATTN  IRVING  L.  KOFSKY 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


01  CY 

ATTN 

P.  FISCHER 

01  CY 

ATTN 

W.F.  CREVIER 

01  CY 

ATTN 

STEVEN  L.  GUTSCHE 

01  CY 

ATTN 

R.  BOGUSCH 

01  CY 

ATTN 

R.  HENDRICK 

01  CY 

ATTN 

RALPH  KILB 

01  CY 

ATTN 

DAVE  SOWLE 

01  CY 

ATTN 

F.  FAJEN 

01  CY 

ATTN 

M.  SCHEIBE 

01  CY 

ATTN 

CONRAD  L.  L0N3MIRE 

01  CY 

ATTN 

B.  WHITE 

01  CY 

ATTN 

R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NEW  MEXICO  87106 
01 CY  R.  STELLINGWERF 
01 CY  M.  ALME 
01 CY  L.  WRIGHT 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01 CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


01  CY 

ATTN 

FORREST  GILMORE 

01  CY 

ATTN 

WILLIAM  B.  WRIGHT 

01  CY 

ATTN 

WILLIAM  J.  KARZAS 

01  CY 

ATTN 

H.  ORY 

01  CY 

ATTN 

C.  MACDONALD 

01  CY 

ATTN 

R.  TURCO 

01  CY 

ATTN 

L.  DeR AND 

01  CY 

ATTN 

W.  TSAI 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01 CY  ATTN  G.  HARDING 
01 CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

UESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01 CY  ATTN  W.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01CV  ATTN  E.C.  FIELD,  JR. 


RANO  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
01 CY  ATTN  ED  BED ROZ IAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01 CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

0 1 C Y  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01CY  ATTN  LEWIS  M.  LINSON 
01CT  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01 CY  ATTN  E . A .  STRAKER 
0 1 C  Y  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GCOuRIOGE  DR. 

MCLEAN,  VA  22102 
01 CY  J.  COCKAYNE 
0 1 C  Y  E.  HYMAN 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 


MENLO 

PARK, 

CA  94025 

01  CY 

ATTN 

J.  CASPER 

01  CY 

ATTN 

DONALD  NEILSON 

01  CY 

ATTN 

ALAN  BURNS 

01  CY 

ATTN 

G.  SMITH 

01  CY 

ATTN 

R.  TSUNODA 

01  CY 

ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHE5NUT 

01  CY 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

01  CY 

ATTN 

J.  VICKREY 

01  CY 

ATTN 

RAY  L.  LEAD A3RAND 

01  CY 

ATTN 

G.  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

01  CY 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V.  GONZALES 

01  CY 

ATTN 

D.  MCDANIEL 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01CY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01 CY  ATTN:  N.  ZABUSKY 


TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01 CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01 CY  ATTN  JOHN  ISE,  JR. 

01 CY  ATTN  JOEL  GAR3ARIN0 

TRW  DEFENSE  8  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  8EACH,  CA  90273 
01CY  ATTN  R.  K.  PLE3UCH 
01CY  ATTN  S.  ALTSCHULER 
0 1 C  Y  ATTN  D.  DEE 
01 C  Y  ATTN  D/  STOCKWELL 
SNTF/1575 


